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The results of the entire Super-Kamiokande-I solar neutrino data are presented. The measured interaction rate 
is 47 ±2% of the rate expected by the standard solar model and 133 ±5% of the rate implied by the SNO charged- 
current interaction rate. There is no evidence for spectral distortion or a time dependent neutrino flux. Together 
with the rates of other experiments, the Super-Kamiokande results imply active solar neutrino oscillations and 
restrict neutrino mixing and mass square difference to lie within the LMA solution area. 



1. Introduction 

Between May 31 st , 1996 and July 15 th , 2001 
Super-Kamiokande (SK) exposed 22,500 tons of 
purified water for 1496 days to solar neutrinos. 
The detector is described elsewhere [Q. Solar 
neutrinos are detected using the Cherenkov light 
emitted by the recoiling electron from neutrino- 
electron elastic scattering. This water Cherenkov 
detection technique offers four advantages: (i) 
the reconstructed direction of the recoiling elec- 
tron is strongly correlated with the neutrino di- 
rection, (ii) the time of each neutrino interaction 
is recorded, (hi) the neutrino spectrum can be 
inferred from the recoil electron spectrum, and 
(iv) neutrino-elastic scattering is sensitive to all 
neutrino flavors, although the cross section for 
and v r is about six to seven times less than for v e . 
However, only high energy solar neutrinos can be 
observed in this fashion; SK analyzes data above 
a threshold of 5 MeV of the total energy of the 
recoil electron. Only extremely few solar neutri- 
nos have an energy that high; they arise from the 
(3 + decay of 8 B ( 8 B neutrinos) and 3 He-proton 
fusion (hep neutrinos) in the sun. The stan- 
dard solar model (SSM) prediction of the flux of 
those rare neutrinos (tpsg = 5.05^' %\ x 10 6 /cm 2 s 
and (j)hep = 9.3 x 10 3 /cm 2 s with no uncertainty 
given [Q) is therefore difficult and beset by large 
uncertainties. 
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Figure 1. Angular Distribution of Solar Neutrino 
Event Candidates. 



2. Solar Neutrino Flux 

SK uses the known position of the sun (at 
the time of a solar neutrino event candidate) 
to statistically separate solar neutrino interac- 
tions from background events. During 1496 
effective days 22,400±800 solar neutrino inter- 
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actions were found in the gray shaded forward 
peak of Figure [|. From Monte Carlo calculations 
based on the SSM, 48,200^;™ are expected. In 
spite of the large neutrino flux uncertainty, the 
elastic scattering rate is significantly less than 
expected from the SSM. The suppression fac- 
tor Data/SSM MC is about 3er less than one: 
0.465±0.005(stat.)t°;°^(syst.)t°; 1 1 ?(SSM). All 
solar neutrinos are born as z/ e in the sun; if all so- 
lar neutrinos detected by SK are i^ e 's as well then a 
v e flux of 2.35±0.02(stat.)±0.08(syst.) x 10 6 /cm 2 s 
leads to the observed interaction rate. 

SNO has measured the charged-current inter- 
action rate of solar v e with deuterium Q (above 
a threshold neutrino energy of about 7 MeV) im- 
plying a v e flux of 1.76t^(stat.)±0.10(syst.) 
(assuming an undistorted 8 B neutrino 
spectrum) and consequently 16,800^'Jqq 
v e interactions in SK. SK has observed 
1.334±0.013(stat.)t°;°g(syst.)j;°;^(SNO) more 
events which is about 4.5<r above one. 

Beyond the 8 B endpoint of » 15 MeV, hep neu- 
trinos dominate the solar neutrino flux. Taking 
into account energy resolution, an optimum hep 
neutrino search window (18 to 21 MeV) was de- 
fined using Monte Carlo. In that window 4.9 ±2.7 
events were observed; one hep neutrino interac- 
tion is expected. This limits the hep flux to be 
less than 73 x 10 3 /cm 2 s (or less than 7.9 times 
the SSM flux) at 90% C.L. 

3. Solar Antineutrinos 

If there are antineutrinos in the solar v flux, 
they could be detected by an inverse (3 reaction 
on hydrogen. Unlike elastic scattering (as shown 
in Figure 0, those events are reconstructed within 
about 6(T of the forward direction), the angu- 
lar distribution of this reaction is essentially flat, 
with a slight backward bias, especially for low en- 
ergy neutrinos. The positron energy is strongly 
correlated with the antineutrino energy. After 
solar neutrino interactions are removed using an 
angular cut of 60°, the dominant background at 
higher energy (above 8 MeV) is due to cosmic 
ray-induced spallation. Figure ^ shows two up- 
per flux limits at 90% C.L.: (i) all events below 
cos #sun < 0.5 are considered to be antineutrino 




Figure 2. Energy-Dependent Upper Limit of 
at 90% C.L. The left panel assumes the 8 B spec- 
trum, the right panel uses monochromatic neu- 
trinos at eight energies. The black data points 
statistically subtract spallation background. 



interactions (ii) spallation background is statisti- 
cally subtracted. The limits of the left plot (in % 
of the SSM 8 B flux) assume antineutrinos with a 
8 B energy spectrum, the limits given in the right 
plot (in /cm 2 s) are for monochromatic antineu- 
trinos of eight different energies. 

4. Time Variations 

Figure |^ shows the time variation of the neu- 
trino flux inferred from the SK rate in bins of 
ten days and 1.5 months. A 7% yearly varia- 
tion is expected due to the 3.5% change in dis- 
tance between sun and earth (assuming the neu- 
trino flux is proportional to the inverse square 
of this distance, that is, the sun is a neutrino 
point source). The data favor such a 7% vari- 
ation over no variation by about 2.5cr. No sig- 
nificant time variation is seen after subtraction 
of this effect. To search for daily variation, the 
data was binned according to solar zenith angle 
9 Z . Figure |I] displays the resulting distribution. 
No significant daily variation is found. The result 
can be summarized by forming the "day/night 
asymmetry" Ajjn = ^d+n) ( w ^h D mean- 
ing "day rate" and N meaning "night rate"): 
A DN = -0.021 ± 0.020(stat.)±g;^(syst.) is con- 
sistent with zero within 0.9<r. 
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Figure 3. Time Variation of the SK Elastic Scattering Flux. The gray data points are measured every 
10 days, the black data points every 1.5 months. The SSM prediction is indicated by the gray band; 
also shown are the SNO charged-current and neutral-current measurements. The black line indicates 
the expected annual 7% flux variation. The right-hand panel combines the 1.5 month bins to search for 
yearly variations; the vertical scale is ten times larger and the zero suppressed. The gray data points 
(open circles) are obtained from the black data points by subtracting the expected 7% variation. 
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Figure 4. Elastic Scattering Flux as a Function 
of Solar Zenith Angle. Note the suppressed zero. 



For precision measurements of the spectrum 
it is necessary to calibrate the absolute energy 
scale (and determine the energy resolution) of 
SK with an outside source. We employed an 
electron linear accelerator to determine the en- 
ergy scale within 0.64% and the resolution within 
2.5% Q. In Figure || the observed recoil elec- 
tron spectrum is compared to the expected spec- 
trum: the ratio Data/SSM effectively plots spec- 
tral distortion. There is no significant spectral 
distortion; the fit to the no-distortion hypothesis 
results in x 2 = 20.2/20 d.o.f. corresponding to 
44.3% C.L. The \ 2 includes three systematic ef- 
fects which will shift all bins of the distribution 
in a correlated way (correlated uncertainties) : the 
uncertainty in (i) the 8 B v spectrum, (ii) the ab- 



4 





., )6 18 . .. 20 

recoil efectron energy in Mev 



Figure 5. Rate Suppression (Upper Panel) and 
Day/Night Asymmetry (Lower Panel) as a Func- 
tion of Energy. The smaller error bars of each 
data point show the statistical, the larger error 
bars the combined uncorrelated uncertainty. The 
gray lines are the ±ler correlated uncertainty due 
to the uncertainty in the 8 B v spectrum and the 
SK energy scale and resolution. All data points 
are dominated by statistical uncertainty except 
the combined rate. 



solute SK energy calibration, and (iii) the SK en- 
ergy resolution. The gray lines of Figure ^ show 
the size of the sum (in quadrature) of (i), (ii) and 
(iii). The same figure also plots the day/night 
asymmetry as a function of recoil electron energy. 
The day/night asymmetry is consistent with zero 
regardless of the recoil electron energy. The plot- 
ted correlated uncertainty is the sum of only (ii) 
and (iii) , since (i) is the same for the day and the 
night. 



lar neutrino experiments |f-|6|]. Solar neutrinos 
are always born as v e ; their conversion proba- 
bility into ia, or v T in vacuum is determined by 
their energy and the distance of flight and the 
oscillation parameters. For the purpose of this 
paper, we assume just two massive neutrinos; in 
that case the oscillation parameters are the mix- 
ing angle 8 (controlling the relationship between 
flavor and mass eigenstates) and the differences 
of the squared masses between the two neutrinos 
Am 2 . 

The presence of matter strongly influences this 
conversion probability. In particular, the mat- 
ter density in the sun is high enough to induce 
resonant conversion [Q at solar neutrino energies 
for a wide range of Am 2 (10~ 8 to 10 _4 eV 2 ). In 
that case, the conversion probability is large or 
even maximal, even if the mixing angle is small. 
The matter density in the earth also affects the 
conversion probability leading to apparent daily 
variations of the solar neutrino flux. 
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6. Effects of Solar Neutrino Oscillations 

Solar neutrino oscillations explain both the ob- 
served SK rate reduction compared to the SSM 
prediction and the observed SK rate enhancement 
relative to the SNO-based expectation (men- 
tioned in section ||) as conversion of v e into other 
active flavors, that is or v r . They also explain 
the rates observed by various radio-chemical so- 



Figure 6. Oscillation Predictions of SK Spec- 
tral Distortion (Left) and Day/Night Asymmetry 
(Right). The spectral distortion is depicted in 
contours of equal x 2 ■ The shown values of both 
the x 2 an d the day/night asymmetry contours in- 
crease logarithmically. Superimposed in gray are 
the usual solution areas. 
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Figure 7. Zenith Spectrum and Four Oscillation Predictions. SMA (dotted line) and VAC (dashed line) 
solutions are disfavored by spectral distortion data, LOW (light gray) solutions by the absence of solar 
zenith angle variation. The LMA solutions (black) fit best. 



Depending on 8 and Am 2 , SK can therefore ob- 
serve distortions of the recoil electron spectrum, 
yearly variations or daily variations. In Figure ^ 
the expected spectral distortion and day/night 
asymmetry are shown for all values of tan 2 9 and 
Am 2 (The spectral distortion is quantified by a 
X 2 method which uses the current SK uncertain- 
ties as input). Overlaid are the traditional os- 
cillation "solutions" , an area of parameter space 
capable of explaining all measured solar neutrino 
rates. The conversion probabilities for the large 
mixing angle (LMA), small mixing angle (SMA) 
and low Am 2 (LOW) solutions are dominated by 
the matter effects mentioned previously, while the 
vacuum solutions (VAC) match the first oscilla- 
tion phase maximum to the solar neutrino energy 
range. LMA, LOW and VAC have a large mixing 



angle close to 7r/4 while the mixing of the SMA 
is small (a few 10" 2 ). SMA and VAC solutions 
predict a distorted spectrum, LMA and LOW so- 
lutions daily variations. 

Figure ^ summarizes the data used for a simul- 
taneous search for spectral distortions and daily 
variations: the zenith angle spectrum has eight 
energy bins and seven solar zenith angle bins be- 
tween 5.5 and 16 MeV. The % 2 for a flat shape hy- 
pothesis is x 2 = 40.5/43 d.o.f. which corresponds 
to 58.1% C.L. Overlaid are the oscillation predic- 
tions (fit to the data by varying freely the total 
8 B and hep neutrino flux) for four selected pa- 
rameter pairs, each representing a solution. The 
strong energy dependence of the flavor conversion 
(inducing spectral distortions) predicted by VAC 
and SMA is not supported by the data. In case 
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Figure 8. SK Excluded (Left) and Allowed (Right) Regions at 95% C.L. The cross hatched areas use 
the SK spectrum, the inside of the dashed lines the shape of the SK day/night spectrum (see Figure ||). 
The shaded areas use the zenith angle spectrum. To get the allowed regions on the right, the shape is 
combined with the SK rate and the SSM constraint on the 8 B flux. 



of LMA and LOW, there is little energy depen- 
dence of the flavor conversion, the apparent spec- 
tral distortion shown in Figure |t] is caused by a 
larger ratio of the hep to the 8 B neutrino flux 
than predicted in the SSM. These two solutions 
can only be distinguished by their daily flux vari- 
ation predictions. LMA shows rapid oscillations 
during the night. These oscillations are washed 
out for the coarse bins of the data. The LOW 
solutions however predict slow oscillations in the 
night. The statistical fluctuations of the observed 
flat distributions would be more consistent with 
slow oscillations of the opposite phase. 



7. Oscillation Analysis 

Figure || shows the strong constraints on so- 
lar neutrino oscillation parameters which are im- 
posed by the SK spectral distortion data (cross 
hatched areas) and the SK zenith angle spec- 
trum (gray areas). A fit to two spectra (day and 
night; see gray dashed lines) partially includes the 
zenith angle variation which is fully included in 
the zenith angle spectrum fit. The analysis is de- 
scribed in ||; the definition of the x 2 employed 
to extract those constraints can be found in ||. 
The areas in the left panel of Figure |^ are ex- 
cluded and use only the (zenith angle) spectrum 
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shape while the right panel shows allowed areas 
based on this shape and the SK rate (and the SSM 
8 B flux). SK rate and spectrum allow only large 
mixing angles. The large daily variation predicted 
between as 10~ 5 and ~ 10 -7 eV 2 splits this "large 
mixing band" into two distinct allowed areas (SK 
LMA and SK LOW/quasi-VAC). The Am 2 lim- 
its of these areas are stronger for the zenith angle 
spectrum compared to the day/night spectrum, 
especially the upper limit for the SK LOW/quasi- 
VAC which is reduced by a factor of two to three. 
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Figure 9. Area Excluded at 95% C.L. by the 
Shape of the SK Zenith Spectrum (Gray Area). 
The exclusion is independent of all SSM v flux 
predictions. Overlaid is the allowed area (inside 
thick black line shaded in light gray) from the 
SNO day/night spectrum M]. The overlap of 
both is shaded dark gray. The cross hatched area 
is allowed by the rate of the Homestake exper- 
iment H and the SSM prediction of that rate. 
Similarly, the gray lines mark the boundary of the 
area allowed by Gallex/GNO @ and SAGE [§. 



Figure 9 compares the constraints of the zenith 
angle spectrum (gray area) with the allowed areas 
of the SNO day/night spectrum ]Io|) , the Homes- 
take experiment and the combined "Gallium ex- 
periments" (Gallex/GNO and SAGE). All SNO 
small angle (tan 2 < 0.1; in particular the one at 
tan 2 — 0.002, which is consistent with the Gal- 
lium experiments), VAC (consistent with Homes- 
take and the Gallium experiments, including the 
one at tan 2 8 — 1), and LOW solutions (con- 
sistent with Homestake and the Gallium exper- 
iments) are excluded. The SNO LMA solutions 
(consistent with all experiments) are restricted 
to a smaller range. The zenith angle spectrum 
also excludes almost the entire Homestake al- 
lowed area (except for LMA and the region just 
below the LOW solutions). 
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Figure 10. Combined Allowed Regions. The light 
gray area on the left is allowed by the shape of the 
SK zenith angle spectrum combined with the SK 
rate and the SNO charged-current rate. The fit 
leading to the dark gray area includes in addition 
the SNO neutral-current rate. The cross hatched 
allowed area results from a combined fit of the SK 
zenith angle spectrum and the Gallium, Homes- 
take and SK rates together (and the SSM v flux 
predictions). The right panel shows the results of 
a fit to all solar neutrino data using the SSM v 
flux predictions except for the 8 B and hep fluxes. 
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8. Combined Oscillation Fits 

Stronger constraints on neutrino oscillation pa- 
rameters can be obtained when the SK results 
are combined with other solar neutrino data. In 
particular, for the first time, a unique solution 
is obtained, when SK results are combined with 
all other solar neutrino data ||. The strongest 
constraints that are independent of SSM neutrino 
flux predictions come from fits which combine the 
SK zenith spectrum shape with the SK and SNO 
rates Q . The left plot of Figure [l0| shows the 
allowed regions of this fit. Superimposed is a com- 
bined fit of the SK zenith spectrum shape and the 
Gallium, Homestake and SK rates which relies on 
the SSM neutrino flux predictions. In either case 
only large mixing angle solutions survive. The 
LMA solutions are favored in either case. When 
all solar neutrino data is combined, the LMA so- 
lutions arc more strongly favored. 

Figure |ll| plots the \ 2 difference A\ 2 as a 
function of only one parameter for various fits. 
(The other parameter is chosen as to minimize 



X 2 -) A\ 2 has the statistics of a \ 2 with one de- 
gree of freedom. SK data by itself strongly dis- 
favors small mixing. LMA solutions are favored, 
however quasi- VAC solutions (between LOW and 
VAC) provide reasonable fits as well with small 
"pockets" of parameter space that fit even bet- 
ter than LMA solutions. The SK-SNO combined 
fit favors the LMA at the 2a level and tightens 
the allowed range in the mixing angle. The SK- 
Gallium-Homestake combined fit disfavors LOW 
and quasi- VAC solutions more strongly, but SMA 
solutions are only disfavored by 2 to 2.5a. The 
all-combined fit favors the LMA by almost 3a 
with small allowed ranges in tan 2 9 and Am 2 . 

9. Conclusion 

Super-Kamiokande solar neutrino precision 
measurements provide the first hint of solar neu- 
trino flavor conversion due to oscillation in a 
unique parameter region, the LMA solution. So- 
lar neutrino oscillation parameters that lead to 
significant distortions of the 8 B neutrino spec- 
trum or solar zenith angle variations are excluded 
by Super-Kamiokande data. 
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Figure 11. Difference in \ 2 as a Function of tan 2 8 
(left) and Am 2 (right). The top plots show only 
the fit to the SK data and the SSM. The bottom 
plots contain the fit to all solar data (black solid 
line), the fit to SK and SNO data (gray solid line) 
and the fit to Homestake, Gallex/GNO, SAGE 
and SK data (and the SSM; gray dotted line). 
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